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Abstract 

In literature, surface tension has been investigated mainly from a thermodynamic standpoint, 
more rarely with kinetic methods. In the present work, surface tension in drops is studied in 
the framework of kinetic theory, starting from the Sutherland approximation to Van Der Waals 
interaction between molecules. Surface tension is calculated as a function of drop radius: it is 
found that it approaches swiftly an asymptotic value, for radii of several times the intermolecular 
distance. This theoretical asymptotic value is compared to experimental values for a few liquids, 
and is found in reasonable agreement. 
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I. INTRODUCTION 



The effect of curvature on surface tension has been investigated widely in the 50's, but 
mainly within the framework of Thermodynamics (see, e.g., jjjEj)- Some studies have taken 
the point of view of kinetic theory 0, Q|. More recent work investigates surface tension, 
in plane geometry, from the point of view of kinetic theory, making use of the Sutherland 
potential |5|" still much work is being done taking the Thermodynamics standpoint (see for 
instance 0, Q, [J). In the present work the problem is studied through kinetic equations 
derived in the following simplifying assumptions: 

1) Free volume method This approximation can also be bettered including the 
molecular exchange with the surrounding atmosphere This effect might form the object 
of subsequent work. 

2) Uncorrelated molecules. This assumption is dictated by two considerations: first, 
the double distribution function (or pair distribution function, as it often referred to) is not 
well known as it can be calculated exactly only for extremely simplified situations. Albeit 
it is almost invariably called upon, still the assumptions made to calculate it are often such 
as to render its benefits very limited; secondly, the importance of correlation can be judged 
from the ratio 7 between the potential energy at the average intermolecular distance and 
the average kinetic energy of molecules [Tol. flll| : 

The intermolecular distance can be estimated, as usual, from the number density n as 
. As will be shown in section IV, where 7 is calculated for several cases, this ratio is of 
order unity, as might have been expected: therefore, correlation between molecules might 
play some role. However, it will be neglected here as a first approximation. 

3) Constant density. The number density will be considered constant with radius: this 
approximation might turn out questionable at the interface, particularly in view of the fact 
that surface tension is, after all, an interface effect. This assumption will be made here, as 
a first approximation, and then reviewed in commenting the results. 

4) Sutherland potential. The intermolecular Van der Waals forces are modeled with 
the Sutherland potential. This assumption could be bettered using the full Lennard- Jones 
potential, however at the cost of greater mathematical difficulties; on the other hand, it can 
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be argued that the error introduced using the Sutherland potential, with suitable parameters, 
is a reasonable price to pay for the mathematical simplification, which allows solution in 
closed form. As will be shown in the last section, the results obtained are in very reasonable 
agreement with experimental data. 

In the following sections, a simple method will be derived within the approximations 
discussed above, to obtain a simple expression for the surface tension as a function of radius. 
As the radius becomes larger, the value of surface tension approaches that of a plane surface. 
Calculations presented for several different liquids, in the large radius limit, are consistent 
with experimental results for plane surfaces. 

II. GOVERNING EQUATIONS 

Consider a spherical droplet of radius R centered in the point O, and a point P at a 
distance r from O, with re [0,R]. A system of spherical coordinates can be defined, with 
origin in P and the O-P direction as the polar depicted in Fig. 1. The (p=0 half-plane 

can be chosen arbitrarily, due to the spherical symmetry of the system. 




Figure 1: Geometry of the problem - section of the droplet on the half-plane tp=0 
Any point Q within the drop that lies on the f=0 half-plane can then be described 
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uniquely with the pair of coordinates p and i? (or p — cos??) as defined in Fig. 1. 

If forces between particles are purely central forces, i.e., depending only on a given power 
a of the distance, then the force due to the attraction of a particle in Q on the particle in 
P can be written as 

F QP = jj (2) 
as it is directed from P to Q. In spherical coordinates, the half circle on the ip=0 half-plane 
depicted in Fig. 1 is described by the equation 



p = yj{prf + (i? 2 - r 2 ) 2 - pr (3) 

or equivalently 

therefore the domain corresponding to the intersection of the drop with the half-plane is 
defined as 



0, \J (pr) 2 + (R 2 — r 2 ) 2 — pr 



/xe [-1,1] ; pe 

or equivalently 

r r (#2_ r 2N_ 2 

P ■:. [0, /.' • /• ; pe 1 



— 1, min \ 1, 



(5) 



(6) 



2rp 

If there are n (ri) particles per unit volume (due to the spherical symmetry the density can 
only depend on the radial position ri, and this can be expressed in terms of r, p and p), in 
a volume element dV (expressed in spherical coordinates) at location (p, p, ip) there are a 
number of particles given by 

ndV = n (ri [r, p, p]) p 2 dpdpd(p (7) 

exerting on the particle in P a force dT 

G 
~P' 

The z-component of dT is 



dT = —ndV = Gn (n) p 2 ~ a dpdpdp (8) 



dT z = dT cos$ = Gn (ri) p 2 ° ' pdpdpdtp (9) 
and in the same vein the component perpendicular to z: 

dT± = dT sin = Gn (n) p 2 ' a y/l - p 2 dpdpdp (10) 
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The first can be integrated readily with respect to ip, yielding 



/ dT z = 2iiGn(r 1 ) p 2 ~ a pdpdp (11) 
Jo 

As for the second, some further discussion is appropriate. To evaluate the surface tension 
one can proceed as follows: imagine cutting the drop in two halves, each one exerting a 
certain amount of attraction on the other. Consider now a molecule that is lying right on 
the edge of one half drop, say at the intersection of the polar axis with the drop surface. 
The component perpendicular to the axis, and lying on on the <p — plane, of the force 
acting on this molecule that is due to an element of volume at location (p,$,ip) is given by 

dJ-±_ cos <p = Gn (ri) p 2 ~ a >Jl — p 2 dpdp cos <pd(p (12) 

The tangential force dFx acting on this molecule due to all elements with coordinates (p,f?), 
is given by the integral over ip of dJ r ±cos(p in one of the halves of the drop, e.g., for 

dF T = / Gn (ri) p 2 ~ a \Jl — p 2 dpdp cos pdip = 2Gn (r x ) p 2 ~ a ^/l — p?dpdp (13) 

To calculate the total tangential force on the molecule, the above expression is to be inte- 
grated over p and p. In the following the density will be considered constant, as discussed 
in point 3) in the introduction. Recalling the expression for the domain in Eq. (0, noted 
that in this case r = R, the following integral is obtained: 

2R r - i r 2R { f-o\ o 




2nG J D p2 ' adp }^ ^^=§7, ^ Q { arCSin V2^ 2R\- \2RJ 

(14) 

Here D is the distance of minimum approach, as will be discussed in the next section. 



III. SOLUTION WITH VAN DER WAALS TYPE FORCES 



Van der Waals interactions are usually represented as a Lennard- Jones potential [1 



This last can be approximated by a Sutherland potential [12j, usually written as 



(15) 



oo r < D 

U s (r) = { (16) 
^ -A(^Y r>D y ' 
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with suitable choice of the (positive) constants uj, A and D. A condition typically imposed 
is to conserve the asymptotic behavior of ULJ, which entails u = 6 and A ■ D 6 = As ■ cr 6 ; 
a second condition can be chosen as follows: the rest position of a particle in the field of 
another particle is the bottom of the potential well, therefore it seems reasonable to choose 
the constants so that this position is the same in the two cases. Therefore, D = a ■ 2 s . With 
the above positions the Sutherland potential becomes 

{oo r < D = a ■ 2e 
(17) 
-ei^Y r>D = o-l\ 

This potential is depicted in Fig. 2. This corresponds to the present model for the force Eq. 
0, with (p=7 and 

G = 6eD 6 (18) 
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Figure 2: Schematized Sutherland potential 

In the Sutherland approximation, D is the distance of minimum approach, the lower limit 
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of the integral in Eq. JHj). Calculating that integral, and recalling the expression for G in 
Eq. (jUj) , the following expression is obtained for the tangential force acting on a molecule: 



3 - I f-D 

Ft = —neD < arcsm , 

2 I \2R 




(19) 



The number of molecule per unit length can be estimated as ns . Therefore the force per 
unit length can be calculated as 

2" 



T sur (R) = ^meD 2 < arcsin 



2R\ 1 \2r) 




(20) 



The graph of the reduced surface tension T r r as a function of ^P, the drop diameter in 
units of D, is shown in Fig. 3 




Figure 3: Reduced surface tension vs. drop diameter in units of D 

As can be gathered from the graph, the surface tension approaches rapidly the limiting 
value of as the drop diameter grows to only a few tens of times D. 
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IV. COMPARISON WITH EXPERIMENTAL DATA 



Limiting values of surface tension for large radii have been calculated from Eq. (|20j) for 
several liquids, and compared with known experimental values for plane surfaces. The results 
are shown in Table 1. As can be gathered from the data reported, the agreement is quite 
good for such a qualitative approach, being always within a factor of 2 to 3 (3.4 for water). 
As the results depend on the parameters of the Sutherland potential, agreement might be 
improved by a different choice, the one choice made here being somewhat arbitrary albeit 
reasonable. In Q. for instance, detailed calculation and comparison with Lennard- Jones 
parameters is reported for several inert gases. 



Table 1: comparison of calculated and experimental tubus of surface tension 




Water, 20°C 
Ethanol, 20°C 



Density 



3.34X10 28 
l.OlxlO 28 



Lennard- Jones | O(r ) 
parameters W j * ~ 3 



2.641 
4.530 



KT 



809.1 ; 
362.6 ! 



1.4 
1.5 



Calculated 
value 
[N/m] 

248xl0 3 
66.5xlO- 3 



Experimental 
value 
[N/m] 

72.75xl0 3 M 
22.75xl0 3 M 



Acetone, 20°C 



0.671X10 28 



4.600 



560.2 i 

L_ 



0.57 



61.5xl0- 3 ! 23.70xl0- 3 M 
62.0xl0- 3 ! 28.88xl0- 3 W 



Benzene, 20°C 



0.677xl(P 8 



5.349 



412.3 ! 2.01 



■+- 



Liquid Oxigen, -183 °C 



2.14X10 28 



3.467 



106.7 | 1.26 



31.2xl0- 3 j 13.2xlO- 3 M 
lS.SxlO" 3 | 6.6xl0 3 CM 



Liquid Nitrogen, -183 °C 1 1.72X10 28 



3.798 



71.4 



0.94 



On the other hand, the density is indeed essentially constant throughout the drop except 

near the boundary, where it is generally found to decrease H Q 0,0, 0, 0, Q- As 

the layers close to the boundary give the main contribution to the integral in Eq. (14), the 
actual density profile should be considered in the integration. This will form the object of 
a subsequent work, for the present the following view will be taken: in the literature, the 
density is always shown to decrease rather abruptly near the interface, taking a relative 
value at the interface of roughly one half the bulk value. Therefore, introducing the actual 
density in the integral will decrease this latter by a factor that is between 0.5 and 1. On 
the other hand, in calculating the number of molecules per unit length at the surface, the 
value of density at the surface must be considered; altogether then, adding the two effects, 
the result changes by a factor of between 0.4 and 0.8, which is the right magnitude to offer 
the correction needed to reconcile calculated values with those found experimentally. 
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V. CONCLUSIONS 



As anticipated in the introduction, aim of this work was to propose a simple equation for 
calculating surface tension in drops as a function of radius. The result is in very reasonable 
agreement with the experiment. This result was obtained in the simplifying assumptions 
discussed in the introduction. Accurate description of the density profile is needed to improve 
confidence in the results yielded by the method proposed. 
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